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In certain two-dimensional time-dependent flows, the braiding of periodic orbits pro- 
vides a way to analyze chaos in the system through application of the Thurston- 
Nielsen classification theorem (TNCT). We expand upon earlier work that introduced 
the application of the TNCT to braiding of almost-cyclic sets, which are individ- 
ual components of almost-invariant sets [Stremler, Ross, Grover, Kumar, Topological 
chaos and periodic braiding of almost-cyclic sets. Physical Review Letters 106 (2011), 
114101]. In this context, almost-cyclic sets are periodic regions in the flow with high 
local residence time that act as stirrers or 'ghost rods' around which the surround- 
ing fluid appears to be stretched and folded. In the present work, we discuss the 
bifurcation of the almost-cyclic sets as a system parameter is varied, which results 
in a sequence of topologically distinct braids. We show that, for Stokes' flow in a 
lid-driven cavity, these various braids give good lower bounds on the topological en- 
tropy over the respective parameter regimes in which they exist. We make the case 
that a topological analysis based on spatiotemporal braiding of almost-cyclic sets can 
be used for analyzing chaos in fluid flows. Hence we further develop a connection 
between set-oriented statistical methods and topological methods, which promises to 
be an important analysis tool in the study of complex systems. 
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When a body of fluid moves, whether it be in the atmosphere, an ocean, or a 
kitchen sink, there are often regions of fluid that move together for an extended 
period of time. As these coherent sets of fluid trace out trajectories in space 
and time, they can be thought of as 'stirring' the surrounding fluid. The en- 
tanglement of these trajectories as they 'braid' around each other is connected 
to the level of chaos present in the fluid system. When the sets return peri- 
odically to their initial positions in a two-dimensional, time-dependent system, 
the entanglement of their trajectories can be used to predict a lower bound on 
the rate of stretching in the surrounding domain. We examine the trajecto- 
ries of such 'almost-cyclic sets' in a lid-driven cavity flow and demonstrate that 
this combination of topological analysis with set-oriented methods can be an 
effective means of predicting chaos. The characterization of the entanglement 
and associated prediction of stretching are achieved through application of the 
Thurston— Nielsen Classification Theorem, which in general classifies the topo- 
logical complexity of homeomorphisms of punctured surfaces. While a rigorous 
lower bound on topological entropy is not available in the absence of exactly 
periodic braiding structures, our approach finds the 'topological skeleton' that 
can be used to get an approximate rate of stretching. 



I. INTRODUCTION 

Qualitative and quantitative analyses of mixing in phase space are often key components 
of understanding the dynamics of a complex system. Our focus is on systems exhibiting 
deterministic behavior, although our discussion is applicable to stochastic systems for which 
the full dynamics can be considered a deterministic 'template' on which an added stochastic 
behavior plays a secondary role. From the perspective of dynamical systems theory, the 
presence of mixing in phase space is synonymous with the existence of chaotic trajectories. 
There has been a significant amount of interest in understanding how such chaotic behavior 
arises, how to detect it, and how to design for or against itP^l 

The field of fluid mechanics has been a particularly productive proving ground for un- 
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derstanding mixing in dynamical systems. Consider, for example, the development and 
growth 4 of the field of chaotic advection, the phenomenon in which advective transport by 
a regular velocity field leads to chaotic particle trajectories. Given an incompressible fluid 
moving with a velocity field V = (u, v, w), one can write the advection equations of motion 
for a passive fluid particle at position x(t) as 

£=V(M). (1) 

Non-autonomous two-dimensional systems or autonomous three-dimensional systems are ca- 
pable of generating chaotic particle trajectories even if the underlying velocity field is not 
chaotic or stochastic. In the dynamical systems vernacular, a mixing system is one that gen- 
erates chaotic trajectories. From the perspective of fluid mechanics, this advective transport 
is typically referred to as 'stirring', with the term 'mixing' usually reserved for homogeniza- 
tion due to the combined effect of stirring and molecular diffusion. It is usually assumed, 
and often correctly so, that dynamical systems mixing is synonymous with fluid mixing, and 
here we will use the terms interchangeably. This approach to understanding and predicting 
mixing has played an important role in the engineering analysis of laminar flows, and it has 
been particularly valuable in the development of microfluidicd^. From the perspective of 
dynamical systems theory, this connection between regular flows and chaotic trajectories was 
previously knowiP. However, it was not until the concept of chaotic advection took root that 
the general scientific community fully realized the ubiquitous existence of relatively simple 
flows capable of generating chaotic transport. 

In a typical analysis of chaotic advection, one requires a detailed spatio-temporal char- 
acterization of transport throughout the flow domain in order to quantify mixing from the 
perspective of ergodic theory using metric quantities. Poincare sections and Lyapunov ex- 
ponents, and the closely-related measure-theoretic or Kolmogorov- Sinai (KS) entropjEEl, are 
generated by accurately tracking individual particles for very long times^l. Quantifications 
of transport using lobe dynamics involve tracking exponentially-stretched interfaces forward 
and backward in time in order to determine stable and unstable manifolds and the area 
(or volume) that they encompas d 10 1 11 1 . Lagrangian Coherent Structures (LCS), which can 
be viewed as generalizations of manifolds, are determined using the finite-time Lyapunov 
exponent field 12 1 13 1 . Each of these methods provides important information about chaos in a 
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domain, but they give information only in those regions of the domain where detailed metric 
calculations can be made, and having incomplete information often prevents analysis. 

In an alternative approach, mixing is quantified using topological, as opposed to metric, 
characteristics of the flow. Applications of topology to analyzing fluid motion date back 
to Helmholt d 14 1 15 1 , but the topological approach we employ here is a modern development^. 
The core idea in this approach is to use the topological classification of a few entangled 
periodic orbits to place a lower bound on the exponential growth rate of material surfaces 
under iteration of the flow. The connection between topology and exponential growth rates 
originated with Adler, Konheim and McAndrew's definition of topological entropy, which 
we will call h, as an analogue of the KS entropy^. There is an enormous body of work 
concerning the dynamical meaning of h and ways to compute and estimate it. Here we focus 
our attention on the main ideas and results that are relevant to the problem at hand. 

The first key idea is the relation between h and the exponential growth rates of various 
quantities induced by iteration of a function. Perhaps the strongest such result is the proof 
by Yomdiif^of Shub's entropy conjecture that h is bounded below by the log of the spectral 
radius of the induced action on homology. That is, the algebraic representation of the 
system topology (i.e., homology) produces a matrix representation of the iterated flow map 
(i.e., the action), and the log of the leading eigenvalue gives a lower bound on h for that 
flow. Newhouse extended this result to the scenario of interest here by showing that for 
C^-diffeomorphisms (i.e., infinitely different iable, invertible mappings) of surfaces, which 
includes all singularity- free deterministic fluid flows, h is equal to the maximal growth rate 
of smooth arcs under the action of the diffeomorphism. He also showed the important result 
that h varies continuously under infinitely smooth perturbations of the action, i.e., in C°°- 
familie^pEU. 

Also of importance here is the related (and earlier) result of BoweipZl that h is bounded 
below by the growth rate of the induced action on the fundamental group, which consists of 
sets of equivalence classes of loops under homotopy, i.e., each set in the fundamental group 
consists of a class of loops that are equivalent under continuous deformation. That is, Bowen 
proved that the growth rate of 'representative loops' gives a lower bound on h. That paper 
also contains what later became a key observation, namely, that by puncturing a surface 
at a periodic orbit and studying the action restricted to the punctured surface, one could 
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get better lower bounds for the entropy using strictly algebraic means. Boweip2 makes the 
prescient remark that Thurston's classification of surface isotopy classes can be used to good 
effect on the punctured surface. 

We thus come to the second key mathematical tool used here: Thurston's theorem on a 
surface punctured at a periodic orbitP^S. Again, this concept has been the subject of a large 
body of work, which we consider here to only a limited extent; See Boylancp2for a excellent 
survey and additional references. The Thurston-Nielsen Classification Theorem (referred to 
hereafter as TNCT) characterizes the isotopy classes of continuous, one-to-one, invertible, 
orientable maps, or homeomorphisms, of punctured surfaces. Since we are discussing the 
application of this work in the context of fluid flows, we will further restrict our discussion to 
differentiable maps of two-dimensional manifolds, i.e. to surface diffeomorphisms. The topo- 
logical complexity in the system is introduced through the relative motions of the punctures 
induced by the flow map. The surfaces under consideration are fc-punctured disks, L^, in 
which the k punctures can correspond to physical obstruction ^ 16 1 281 , periodic orbitd^SEES r 
other stirrer-like objects that entangle the surrounding fluid, such as point vortices^. Since 
it is topology being considered here, the specific shapes of the disk and the punctures do 
not matter. In dynamical systems applications, these surfaces are typically two-dimensional 
physical space and the maps act in time, but under certain conditions maps can also be 
constructed in three-dimensional spac^l. 

An isotopy class consists of all the topologically equivalent diffeomorphisms that can be 
mapped to each other through a 'fixed stirrer' diffeomorphism in which the punctures are 
fixed relative to each other during the mapping. According to the TNCT, each isotopy class 
contains a representative homeomorphism the Thurston Nielsen (TN) representative, 
that is one of three possible types: finite order (FO), pseudo-Anosov (PA), or reducible. If 
$ is FO, then its n th iterate is the identity map in which the punctures are held fixed and all 
points in map to themselves, so no net stretching occurs. All fixed stirrer diffeomorphisms 
are isotopic to the identity. If $ is PA, then it stretches and contracts everywhere on the 
surface by factors Atn > 1 and A^. Iteration of a PA $ results in exponential stretching 
of smooth arcs with corresponding topological entropy /itn = In Atn- There are several 
methods by which one can compute Atn; for our analysis, we make use of the open-source 
code written by HalPS, which implements the Bestvina-Handel algorithmP^. A necessary 
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(but not sufficient) condition that an isotopy class contain a PA $ is the existence of at least 
three punctures in the surface. Finally, for the reducible case, $ leaves a family of curves 
invariant, and these curves delimit regions in which $ is either FO or PA. Each isotopy class 
contains only one TN representative, so the isotopy class itself can be referred to as being of 
FO, PA, or reducible type. We focus here on cases for which the isotopy class is of PA type. 

The power of Thurston's theorem in the analysis of dynamical systems comes from the 
connection between the TN representative $ and the complexity of a flow / within the same 
isotopy class. By Handel's isotopy stability theorem^, if $ is PA and / is isotopic to $, then 
there exists a compact, /-invariant set Y C and a continuous, onto mapping a : Y — >■ 
such that af = <§>a. Thus, the complex dynamics of a PA $ are preserved by / in (some 
subset of) the domain. The map a may be many-to-one, so that the dynamics of / may 
be more complicated than those of $, but they are never less complicated. Therefore, the 
existence of a PA $ places a lower bound on the complexity of any flow map / from the 
same isotopy class. Non-trivial material lines, such as curves that encircle two punctures or 
that connect punctures with each other and/or the disk boundary, will grow exponentially 
in length under n iterates of / according to L ~ A/ n , with associated topological entropy 
hf = In A/, and the TNCT establishes that Xf > Atn (and hf > /&tn)« Unfortunately, there 
are no restrictions on the size of Y , so the complicated dynamics may exist on a subset with 
Lebesgue measure zero. However, the complexity associated with the exponential stretching 
and folding of material lines cannot be removed by continuous perturbation of the domain 
that maintains the orbit topology. This topological chaos is thus "built in" to the system 
due to the orbit topology of only a finite number of punctures. 

The initial application of the TNCT to fluid mixing examined a system of three cylindrical 
rods moving slowly on periodic trajectories through a viscous fluid in a cylindrical domaiif^. 
The three rods exchanged positions pair-wise on circular trajectories. Two stirring protocols 
were considered: one motion of FO type and one of PA type. In the FO protocol, all mo- 
tions involved clockwise interchanges. In the PA protocol, the exchanges alternated between 
clockwise and counterclockwise, producing figure-eight rod trajectories. Although these two 
motions were energetically equivalent (in the Stokes' flow limit), the mixing produced by 
the PA motion was clearly better than that produced by the FO motion. Computational 
analysis of this flow system 28 confirmed that the lower bound predicted by the TNCT gives 
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an excellent representation of the actual stretching rate caused by the PA motion. Fur- 
thermore, a substantial subset of the domain exhibited exponential stretching and folding, 
demonstrating that the dynamics of the TN representative can indeed be relevant to the 
analysis of a realistic fluid system. 

Artin's braid groupP^ provides a useful framework for representing the trajectories of 
punctures under the action of the flow. Each of the puncture trajectories is identified with 
a strand in a "physical braid" . The temporal reordering of these strands dictates the topo- 
logical structureof the physical braid. The generator Oi represents a clockwise exchange of 
strands i and i + 1, and (our shorthand notation for a^ 1 ) represents a counterclockwise 
exchange. A braid can be identified via a "braid word" of generator "letters", which we 
read from right to left as time progresses. For example, the braid on three strands that 
we discuss in Section [II] is represented by the braid word cr_i<7 2 . Since the braid word, or 
the "mathematical braid", encodes only the topology of the braid crossings, each isotopy 
class on a punctured orientable two-dimensional surface is associated with one such braid. 
Thus, we can refer to a braid as being of FO, PA, or reducible type. A brief overview of 
the connection between braid theory and topological chaos is given in Section II A , and a 
general introduction to braid theory can be found in, e.g., | 37 | 38 [ 

From the mathematical point of view it is clear that the topological analysis of a fluid 
system can be based on the trajectories of periodic orbits, even when these orbits are not asso- 
ciated directly with stirring rods. Such orbits can still appear to be 'stirring' the surrounding 
fluid, and they have been thus been termed ghost rods 29 . In many cases, consideration of 
ghost rods is essential to understanding the presence of topological chaos in a system. For ex- 
ample, a single physical rod moving on an epicyclic trajectory in a two-dimensional domaiiP^ 
produces a braid on one strand, which is in an isotopy class of FO type. However, it is clear 
from direct observation that stirring a Stokes' flow with this motion produces significant 
stretching and folding of material lines. This complexity can be explained using the TNCT 
by noting that there are two important fixed points, or ghost rods, in the domain about 
which the moving rod winds. The braid on three strands corresponding to the space-time 
trajectories of the two fixed ghost rods and the moving rod is in an isotopy class of PA-type, 
and thus predicts the presence of chaos in the domain. The trajectories of moving ghost rods 
can also produce topological chaos, even when there are no physical rods in the system, such 
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as in the lid-driven cavity flow discussed in Section [II] that forms the basis for our present 
discussion. 

Identifying braiding of ghost rods is an important step in applying the concepts of topolog- 
ical chaos to a broad range of fluid systems with no physical stirring rods. However, finding 
appropriate periodic orbits is a challenging task. One way of getting around this difficulty 
is to calculate the long-time braids generated by random collections of initial point^^. 
The topological entropy is estimated from these computed braid generators. This approach 
removes the need to identify periodic orbits in the flow but, since the braids being considered 
here are aperiodic, the TNCT cannot be invoked in the prediction of a lower bound. Instead, 
relatively long-time calculations are needed. It has also been found that a random selection 
of trajectories in the flow typically leads to a poor estimate of overall system behavior™ so 
a large statistical representation is required in this approach. 

We have recently considered a different approach to examining topological chaos in flows 
without physical rods and with no discernible low-order periodic orbits on which to base a 
topological analysis 41 . In this view, the domain is divided into distinct subsets such that 
there is a very small probability that typical trajectories beginning in each subset will leave 
this subset in a short time. These Almost- Invariant Sets (AlSsjP' can be determined from 
the eigenspectrum of the discretized Perron-Frobenius transfer operator via a set-oriented 



approach, which we review in Section HI A In some cases, disconnected components of an 
AIS correspond to almost-periodic regions, or what has been referred to as Almost- Cyclic 
Sets (ACSs) 42 . Since the ACSs consist of trajectories that move together for a relatively long 
time, the region of the domain corresponding to an ACS can be identified as a (leaky) ghost 
rod, and a representative trajectory from each ACS can be used to reveal the underlying 
braid structure. Furthermore, since these ACSs contain almost-periodic orbits, the aperiodic 
braid they generate can be approximated by periodic continuation to a time-periodic braid. 
This time-periodic braid can be characterized using the TNCT, and an estimate of the 
topological entropy can hence be determined. The topological entropy given by the TNCT 
is no longer a strict lower bound, since the time-periodic braid is only an approximation of 
the true dynamical structure. However, it was demonstrated^ that this estimate can give a 
good representation of the flow. 

Our earlier work^H considers creeping flow in a lid-driven cavity and starts with parameters 



for which there exist three periodic orbits, which we refer to as the reference case. For small 
perturbations from the reference case, consideration of ACS braiding gives a clear explanation 
for the topological entropy in this flow. However, for certain larger perturbations from the 
reference case, we observe a drop in the numerically computed topological entropy below 
the lower bound predicted for the reference case. This paper extends the previous analysis 
to consider these larger perturbations, which leads us to our main observation. As we vary 
away from the reference case, the AIS/ACS structure appears to bifurcate, leading to a 
sequence of topologically distinct braids with differing numbers of strands. The motion in 
Dk for each of these new braids on k strands belongs to an isotopy class of PA type, and the 
corresponding value of Atn gives a correct lower bound estimate on the topological entropy 
for the corresponding parameter value. Hence, we give further evidence that these almost- 
cyclic sets are natural objects on which to base an application of the TNCT. We assert 
that this work marks an important step in making the ghost rod methodology applicable 
to realistic fluid flows with arbitrary time-dependence when low-order periodic, braiding 
orbits are difficult to identify or do not exist. We also conjecture that this generalization is 
applicable to a wide variety of dynamical systems, not just fluid systems. 

II. BRAIDING OF PERIODIC ORBITS IN A LID-DRIVEN CAVITY 
FLOW 

A. The reference case 

The fluid system model that we analyze in this work is a simplified version of a systemP^ 
that allows for exact solutions. This relatively simple fluid system has an easily- visualized 
piecewise-steady velocity field, but due to the imposed time dependence it exhibits compli- 
cated dynamics. 

We consider a two-dimensional lid-driven cavity flow in an infinitely-wide cavity with 
height 2b. Under the assumption of Stokes' flow, the stream function ip(x,y) defined by 




(2) 



satisfies the two-dimensional biharmonic equation, 



V 2 V 2 tp(x,y) = 0. 



(3) 
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We assume that the flow is driven by prescribed tangential velocities on the top and bottom 
boundaries at y — ±b. The piecewise steady tangential velocities we take to be 



u{x ) b) = —u(x, —b) = 



Ui sin O +C/ 2sin(^) when m s < t < (re + 1)t//2 

-U x sin (*f ) + U 2 sin (^) when (n + l)r//2 < t < (n + l)r f 

(4) 

for integer n, where T/ is the time period of the system. The solution to Eq. ^ subject to 
the boundary conditions in Eq. Q is 



ip(x,y,t) 



UiCMy) sin (f ) + U 2 C 2 f 2 (y) sin ( 2 -f) nr f <t<(n + l)r f /2 

-t/idAO/) sin (f ) + £/ 2 C 2 / 2 (y) sin (?f ) (re + 1)^/2 < t < (re + 1)^, 

(5a) 

where 

//c(?/) cosh sum sum cosh fob) 

a \ a J \ a J a \ a J \ a J 



and 



Ck 



a 2 \ a f 2knb\ 



2k7T 2 b 



sinh + 1 



(5c) 



_2k7rb 

The spatial symmetry in the boundary conditions produces a vertical streamline in the flow 
at x = na, and without any loss of generality we can restrict our attention to the bounded 
two-dimensional rectangular domain 

D = {(x, y) : < x < a, -b < y < b}. (6) 

By the time periodicity of the boundary conditions, the flow pattern is reflected about the 
line x = a/2 every t//2 units of time. 

In Section 3 we will often refer to the one-parameter family of stroboscopic global Poincare 
maps, 

4>l +Tf :D^D, (7) 

where (j)^ Ts is the solution diffeomorphism (i.e., flow map) from time t to t + Tf\ that is, 
(^ +r/ represents the motion of all fluid particles over one period of the velocity field starting 
at time t G [0,t/]. Here we treat the initial time t as a parameter, which we can also think 
of as an initial phase (belonging to S 1 ). 
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We will focus on analyzing behavior in a domain with aspect ratio a /(2b) = 3, such 
as is shown in figure [T] The structure of the flow for the reference case depends on the 
value of the parameters a, C/i, L^, and r/. The ratio U2/U1 dictates the streamline pattern, 
and the relationship between U\ (say) and 77 determines how much time a particle spends 
traveling along a streamline before the flow is 'blinked' to another streamline pattern. For 
a domain aspect ratio of aj (26) = 3, we have the following properties when we take a = 6, 
U 2 /U x « 0.841, [/1 « 0.1055, and 77 = 1: 

1) There exist three points x L , and x#, such that is at the center of the domain 
(x, y) = (a/2, 0), and xl and x# are located symmetrically about x = a/2 along y = 0. 

2) For < t < T//2, x L is a fixed point, and and xr exchange their positions while 
moving clockwise along their shared streamline. 

3) Since the flow pattern is reflected about x = a/2 at t = T//2, x^ is a fixed point for 
Tf/2<t<Tf, and xc and x L exchange their positions during this time period while 
moving counterclockwise. 

After three periods of the flow, the points x^, x^ and x# return to their original positions. 
We refer to this choice of parameters as the reference case, for which the trajectories of these 
three period-3 points form the (2+l)-dimensional space-time braid on three strands, a_i<7 2 , 
shown in figureflVc). We determine this braid structure by projecting the trajectory crossings 




(a) (b) (c) 



FIG. 1. (a-b) Representative streamlines given by Eq. ^ for the reference case with a/ (26) = 3 
and U2/U1 ~ 0.841 for (a) < t < Tf/2 and (b) Tf/2 < t < Tf. Points marked by solid circles are 
stagnation points in the flow. If a = 6 and U\ « 0.1055, the fluid particles starting at the open 
circles have exchanged positions after At = Tf/2 = 1/2. (c) The braid diagram generated by the 
periodic orbit trajectories in the reference case. 
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onto the x-axis (which gives the 'physical braid representation') and follow the braid group 
labeling convention discussed in Section 1. 

In the reference case, the determination of Atn can be described via the Burau representation!^. 
In this representation, a braid on TV strands is identified with a (N — 1) x (N — 1) matrix. 
Each of the braid generators has a matrix representation, and the braid matrix is given by 
the product of the generator matrices. The entries of Burau matrices are Laurent polynomi- 
als in some variable i, which is taken to be —1 here. The Burau matrices can be heuristically 
understood as Markov matrices that track the evolution of material lines connecting the 
points. For instance, in figure [2| the generator a 2 stretches the line I to approximately look 
like the sum of lines I and II while leaving the line II intact, giving the relations I'=I+II and 
IF =11. The Burau matrix corresponding to a 2 is thus 



M 2 = 



1 
1 1 



(8) 



Similarly, for <r_i the Burau matrix is, 



M_i 




(9) 



The Burau matrix for the braid a-i(J2 is then M = M_i Mi- According to the TNCT, when 
N = 3 the stretching rate produced by the TN representative $ is the same as in the linear 
map represented by the Burau matrix. That is, for the reference case Atn is given by the 
dominant eigenvalue of M, so that Atn = (3 + v / 5)/2. In order to distinguish between the 



(a) O 



2 

O 



(b) O-j— O 
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-o 
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-o ■=> 6— o 
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FIG. 2. The stretching of material lines under the action of the two generators in the reference 
case: (a) the action of 02, (b) the action of cr_i, and (c) the actual stretching of a material line 
that is initially along the x axis in the reference case. 
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stretching produced by topologically distinct flows with TV punctures, we will refer to the 
entropy in this reference case as 

^TN,3 = log(A T N) ~ 0.962. 
For TV > 3, the dominant eigenvalue of the Burau matrix only provides a lower bound for 
Atn, and other methods must be used instead, such as the train-tracks algorithm ! 33 1 34 1 or the 
encoding of loop d 4Q|45 J . 

The actual topological entropy produced by the flow map, hf, can be determined by com- 
puting the asymptotic stretching rate of topologically non-trivial lined^, such as lines that 
join a periodic point with the outer boundary. We estimate hf by following two orthogonal 
lines that are initially along the lines x = a/2 and y — 0, respectively, as they are stretched 
over 6~10 periods of the flow. For numerical reasons these lines are not extended all the way 
to the boundary, but instead are terminated a distance e < a/200 from the boundary. For 
this reference case, these computations give hf « 0.97, which is reasonably well represented 
by the lower bound, /itn,3- 



B. Perturbation of the reference case 

For the analysis discussed here, we keep the ratio U2/U1 fixed while changing the value 
of Tf away from unity. This perturbation increases (for 77 > 1) or decreases (for 77 < 1) the 
amount of rotation in the domain before switching, which has a direct effect on the stretching 
of material lines. The computed topological entropy for the flow is shown in figure [3] for a 
range of 77 values. The entropy varies smoothly, as expected from the results of Newhous^. 
The overall trend in hf is logical, as increasing Tf increases the energy added to the flow 
during each period. However, the rate of increase in hf is not uniform, which we explain by 
considering the topology of the flow map. 



are 



In the reference case with Tf = 1, the three period-3 points discussed in section II 
parabolic points, which are structurally unstable. Under small perturbations for Tf > 1, each 
parabolic point bifurcates into a pair of periodic saddles and periodic centers, as illustrated in 
figure[4| The stable and unstable manifolds that pass through the saddles (i.e., the hyperbolic 
points) form barriers to transport, so that fluid within these regions moves together for a 
number of periods. These 'ghost rod' structures are 'leaky', in that there is a small amount 
of transport across their boundaries due to the transverse intersections of the stable and 
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0.85 0.90 0.95 Tf 1.00 1.05 



FIG. 3. Variation of the actual topological entropy of the fluid, hf, as a function of Tf. Also shown is the 
lower bound /itn,3 for C-i&2: the braid on 3 strands corresponding to the reference point with Tf = 1. 



unstable manifolds, form lobes. The amount of leakage per period is the amount of phase 
space enclosed by the lob^R In figure |I](b,c), we show the manifold structure for the center 
'leaky ghost rod' for Tf = 1.01 (panel b, where the the transversal intersection is present, 
but not evident) and Tf = 1.04 (panel c). The four period-3 points (two hyperbolic and two 
elliptic) that form the structure of a single 'leaky ghost rod' produce a braid on 12 strands. 
The four strands corresponding to each individual ghost rod structure simply twist around 
each other, and as a result the braid on 12 strands is reducible to a braid on 3 strands. This 
mathematical braid on 3 strands is identical to the mathematical braid from the reference 
case. Thus, the value of /itn ~ 0.962 predicted by the TNCT for the reference case remains 
the lower bound on hf for the perturbed flow we have considered here with Tf > 1. Clearly 
there is additional stretching in these cases that is not captured by the TN representative, 
which suggests the presence of additional ghost rods. We do not explore this increased 
complexity here. 

When the value of Tf is decreased below 1, the parabolic points disappear and no low- 
order periodic points are found in the flow. We have searched numerically for periodic orbits 
by iteratively mapping regions of the flow forward and backward in time. For example, 
the Poincare section shown in figure [4](d) for Tf ~ 0.99 appears to be completely chaotic. 
Thus, there are no apparent periodic 'ghost rods' available to form a braid and predict a 
lower bound on the topological entropy of the flow. However, figure [3] shows that Htn,3 is a 
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FIG. 4. (a) Poincare section from the map (/>Q f for rj = 1.01. Note the 'leaky ghost rods' generated by 
perturbation of the parabolic periodic points from the reference case, (b, c) Close-ups of the domain near 
(x, y) = (a/2, 0), which show the structure of the center 'ghost rod' consisting of two saddle and two centers 
for (b) Tf — 1.01 and (c) Tf = 1.04. In (b) the stable and unstable manifolds nearly overlap, while in (c) 
they clearly intersect transversally, leading to significant 'leaking' in and out of the region bounded by them, 
(d) Poincare section from the map (j) T Q f for Tf « 0.99. The phase space appears featureless and is devoid of 
low-order periodic orbits that form non-trivial braids. 



lower bound for hf for approximately a 7% perturbation in 77™ In the discussion below, we 
consider the validity of the topological entropy predicted by TNCT even if we cannot identify 
exactly periodic orbits that produce the expected braiding motion. We accomplish this by 
using a transfer operator approach to reveal phase space structures that braid non-trivially 
and persist under perturbations in the 77 values. 
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III. ALMOST-INVARIANT AND ALMOST-CYCLIC SETS 



A. Computation of almost- invariant and almost-cyclic sets 

We use a set-oriented method to compute AISs for the lid-driven cavity flow system. Our 
aim is to partition the phase space into a given number (say k) of sets {^4i, A 2) ....A^} such 
that the phase space transport between these sets is very unlikely We also want these sets to 
be important statistically with respect to the long term dynamics of the system. Following^, 
this can be formulated as follows. We define the invariance of a set Ai as 

p ^ Ai) = KA) • (10) 

This quantity represents the probability (according to an invariant measure /i) of a point 

in set Ai being mapped back to Ai over one iteration of a map / = (pl +Tf . We want to 

k 

maximize the quantity ^^p^Aj), with the constraint that fi(Ai) is not too small compared 

i=i 

to 1. Without loss of generality we can restrict to partitions {Ai, A 2 , ....A/,} such that each 
member Aj of is a union of sets in {£>i, B 2) . . . , B n } ) i.e., each Aj = U^/S^, for some set of 
box indices 7c{l,2,....,n}. Using the Ulam-Galerkin methocP^^, we define the transition 
matrix P tjTf with entries 

_ m(f-\B i )nB ] ) 
j ~ rn{Bj) ' 1 } 

where Si, . . . , B n are the boxes in the covering and m is the normalized Lebesgue measure, 

which coincides with the phase space volume measure (see figure [5]). In our computations, 

all boxes will have the same measure, i.e., m(Bi) = m(Bj) for all z, j. 

This stochastic matrix P tjTf is a discretized approximation of the Perron-Frobenius oper- 
ator of the map / = (j)^ Tf that may be viewed as a transition matrix of an n-state Markov 
chain 51 . The first left eigenvector v\ of Pt, T/ , corresponding to an eigenvalue Ai = 1, is the 
invariant distribution of the system, i.e., the discretized version of the invariant measure /jpl. 

Following^, we form a new reversible matrix i?t, T/ , which is a stochastic matrix having 
only real eigenvalue^ that satisfies important properties related to almost-invariancd^. 
First, we define the reverse time transition matrix Pt, Tf i whose elements are given by 



Pij — Pji 
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FIG. 5. Computation of the transition matrix P tjT by a set-oriented method. The phase space is divided 
into boxes, {Bi, B 2l . . . , B n }. Box Pj at the final time t + t/ is mapped (backwards) to f~ 1 (Bj) at the 

t + Tf 

initial time t, where f = 4> t . The value of the entry p^- is the fraction of box P^ that is mapped into box 
by /. 



Here, by ^ we mean the zth component of the vector v\, i.e., the value of v\ assigned to the 
box Bi. (For additional accuracy, one can directly calculate P t ^ Tf as P t? _ r/ 47 .) The reversible 
matrix R tTf is 

We use the left eigenvectors corresponding to eigenvalues of i?t, T/ , to detect AISs. If 
an AIS (over one period-77 map, <^ +Tf ) can be further decomposed into period- TV subsets, 
then those subsets are ACSs of the flow of period iVY/, and are hence AISs under the map 

{<t>\ +TI ) N = <t>\ +NT} . 

For this system, we approximate the dynamics on D by covering it with a collection of 
n = 19200 equally sized square boxes (aspect ratio 1), i.e., 240 boxes in the x direction 
times 80 boxes in the y direction. We take 100 uniformly distributed points in each box (a 
10x10 grid), perform a forward iteration for each point for one period of the flow 77, and 
monitor the boxes reached by the iteration. The transition probability from a source box to 
a destination box is measured by the ratio of the number of points from the source box that 
reach the destination box in one iteration step. 

All subsequent computations were verified to depend weakly on the number of points 
per box (for values greater than 100) by comparing the results with those obtained when 
using 2025 points per box (a 45x45 grid) instead of 100, with the difference in the various 
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eigenvalues for the two cases being less than 0.1% for the range of parameters considered here. 
Similarly, a finer discretization of phase space, obtained by using 30000 boxes (300 x 100 
boxes) was used to verify the convergence of eigenvalues. It was also verified that finer 
discretization does not change the order of the eigenvectors. 

The system being considered here is a canonical Hamiltonian system, and thus the in- 
variant measure, /x, is the same as volume measure, m. As our boxes are of equal sizes, 
the entries of v± should be approximately equal. We note that the invariant measure is not 
unique in our case, since a point measure at the periodic points (if any) will also be invariant. 
However, due to the discretization procedure, we recover the unique absolutely continuous 
measure. 

We use the sign structure of the Vk, k > 2, to detect AlS d 47 1 54 1 . For instance, for a given 
scalar value c (chosen as described below) , the set given by the union of regions in phase space 
such that the component of v 2 is greater than c corresponds to one AIS, while the union of 
regions such that the component of v 2 is less than c corresponds to its complement, another 
AIS. More precisely, the two AISs are defined as follows. Let l\ — {% : v 2 ^ > c, 1 < i < n} 
and I 2 = {i : v 2 ,i < c, 1 < i < n} ) then A\ = U^i?^ and A 2 = Ui e i 2 Bi. The value of 
c is chosen so as to maximize min^p^Ai), p^(A 2 )). For the majority of cases, we will take 
c = 0, but to maintain consistency across figures, for non-zero c we plot (v 2 ^ — c) values 
while identifying the AISs, so that the zero contour on the figure always refers to boundaries 
delineating the AISs. 

It has been showii^ that if p^A^ < p M (A n ), then p(A m ) < p(A n ). It has also been 
shown 52 that for any A m) such that p{A m ) < |, there exist the following bounds on transport, 

pMm) e (l - \/2(l-A 2 ), |(1 + A 2 )) . (12) 

Before moving to an investigation with tj < 1, it is instructive to study the transport 
in our system for tj > 1, since we can compare the set-oriented method results with those 
obtained from lobe dynamics. In figure [HJ we show the stable and unstable manifolds of the 
period-3 saddle points for the case t/ « 1.06. On the right, we show the AIS decomposition 
into two sets (red and blue) near the point (x,y) = (3,0). We are interested in finding the 
transport out of the red region under / 3 , where / = (f>l +Tf . For the transition matrix for 
/, A 2 = 0.9967, and thus for f 3 we have A 2 = (0.9967) 3 = 0.990. Hence the lower and 
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FIG. 6. (Left) The stable (magenta) and unstable (green) manifolds of the period-3 saddle points for 
Tf « 1.06. (Right) The manifolds superimposed with a two-set decomposition of the phase space, shown 
near the point (x,y) — (3,0). 



upper bounds on invariance (12) tell us that, under /, p^{Ai) G (0.857,0.995), where A\ 



and A 2 are defined as before, and fi(Ai) < ~. Lobe dynamics theor^^^S tells us that the 
amount of phase space transported out of a separatrix-bounded region (boundaries given by 
stable and unstable manifolds from saddle to primary intersection points) is given by the 
size of the lobe. Using box counting, we find that one lobe consists of approximately 56 
boxes. There are 2 lobes, one on each side, and hence the approximate transport out of 
the separatrix-bounded region under f 3 equals 112 boxes. Dividing this area by the size of 
the separatrix-bounded region (=1158 boxes), we find that the invariance predicted by lobe 
dynamics is = 1 — 112/1158 ~ 0.903, well within the bounds provided by A 2 . 

In figure [7J we show the AISs based on v 2 for 77 = 0.99. The collection of the three 







FIG. 7. AIS structure for rj « 0.99 based on v 2 . The zero contour (black) is the boundary between the 
two almost invariant sets. Compare with the Poincare section for the same Tf shown in figure Qd) 



prominent regions corresponding to the set A\ (i.e., the union of regions where v 2 > 0) forms 
an AIS. This almost invariant set consists of three ACS components, which form a 3-stranded 
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braid, as discussed next. The eigenvector that has TV ACS components, and hence forms 
an TV-stranded braid, is referred to as an TV-stranded eigenvector in our discussion. Each 
(positive) lower eigenvalue has an associated eigenvector whose zero contour isolates other 
AIS, which are more 'leaky', i.e., the lower the eigenvalue, the lower is the invariance of the 
associated AIS. For reference, we show the next 4 eigenvectors in figure [8j which tend to 
reveal smaller scale structures. 




(a)v 3 



(b> 4 




(d)v 6 



FIG. 8. The eigenvectors 3 through 6 of the reversibilized transfer operator, Ro, Tf , f° r r f ~ 0.99. Here red 
denotes highly positive regions, blue denotes highly negative regions, and green/yellow are close to 0. 



B. Braiding of almost-cyclic sets 



Recall from Section |IIB| that while the three period-3 fixed points cease to exist for 
Tf < 1 (see figure [l](d)), we can see from figure [7] that the three ACS, each of period-3, 
still exist in the same region in phase space for r/ « 0.99. In what follows, we examine the 
trajectories of these sets in space-time in order to determine if the corresponding braid is of 
PA type. We use Tf w 0.975 for our analysis in this case. To construct the physical braid 
for the ACS trajectories, we form a discretized family of time-shifted versions of the transfer 
operator Pt ilTf for tj = (i — l)r//fc, where z = l,2,...,fc. That is, we are considering transfer 
operators corresponding to a family of period- 77 flow maps, ([7]), i.e., 



>u .D->D, 



fti+ T f . 
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of different initial phases U G [0, T/). 

For each of the corresponding i?t-, T/ , we find the almost-invariant structure based on 
the second eigenvector as mentioned previously, which reveals the three period-3 ACS for 
Tf w 0.975. We take = 20, and in figure [9](a-f) we show these structures for the first 
half-period of the flow. This physical braid for this half-period of motion corresponds to 
the mathematical braid generator a 2 . Similarly, for the second half of the time-period of the 
flow (not shown), the stirring protocol is the one given by the braid generator 0-\. Hence, 
the physical braid formed by the ACS trajectories for tj ~ 0.975 is represented by the 
mathematical braid cr_i<T2, which is identical to the mathematical braid formed by periodic 
points for r/ = 1 or by elliptical islands for r/ > 1. Three periods of the flow (one period of 
the braid) are shown in (2+l)-dimensional space-time in figure [9](g), which can be seen to 
be isotopic to the braid shown in figure [jjc). Hence the lower bound on topological entropy 
is again given by /itn ~ 0.962, and we can see from figure [3] that this lower bound is still 
valid. 






FIG. 9. Braiding in (2+l)-dimensional space-time via three almost-cyclic sets for Tf « 0.975, a parameter 
value for which corresponding periodic points no longer exist (enhanced online). The 2nd eigenvector, v^, of 
Rt, Tf is shown for different phases for the first half-period, specifically t/rf = (a) 0, (b) 0.1, (c) 0.2, (d) 0.3, 
(e) 0.4, (f) 0.5. (g) The 3 ACSs shown braiding in space-time for 3 periods of flow (enhanced online). The 
braid is isotopic to the braid shown in figure [TJc). 
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C. Persistence and bifurcation of almost- invariant sets 



As the parameter Tf is varied, one expects to see continuous variation in the eigenvalues 
of i?£, rr Since the eigenvalues of R t , Tf are independent of the phase £, we can consider the 



initial phase case, Ro^ Tf . In figure 10, we plot the first few eigenvalues of i2o jT/ in dependence 
on the parameter 77, similar tcP[ Although several 'crossings' of eigenvalue curves seem 




f j) 




FIG. 10. Almost-invariant set structure via continuation of the 3-strand eigenvector found at (a) for various 
values of rj. The structures shown in (a) and (b) can clearly be seen to again have 3 strands. The structure 
in (c) consists of 16 strands, while (d) can be seen to have 13 strands, (e) has 10 strands while (f) consists 
of 8 clearly identifiable strands. Also shown are the four largest eigenvalues of the reversibilized discretized 
transfer operator Ro, Tf (except 1), A2 > ••• > A5, which are colored according to their ordering. One of 
these is singled out as the A* branch (shown with large squares), which is the eigenvalue corresponding to 
vl f , the continuous eigenvector family of interest (see text for details) . 

visible, care must be taken to determine if these are genuine crossings or simply close ap- 
proaches, as eigenvalues generically 'avoid' crossings if there is no symmetry present^. Our 
concern is not necessarily to resolve families of eigenvalue curves, i.e., to determine whether 
genuine eigenvalue curve crossings occur. Instead, we are interested in determining families 
of eigenvectors by the method of continuation. In particular, we want to determine the 
family of eigenvectors that generate a consistent trend of PA braiding for tj < 1. 

We find that such a structure can be found by continuity of the branch of eigenvectors 
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of i?o,7Y that are connected to v 2 of i? ,i, since v 2 of i? ,i captures the braiding structure 
correctly for the reference case. The continuation of this 3-strand eigenvector for t/ < 1 is 
carried out as follows. 

First, we discretize the parameter space 0.8 < t/ < 1 into intervals of size At/ ~ 0.0021. 
Then we calculate the reversibilized Perron- Frobenius operator Ro^ Tf for each of those param- 
eter values. We compute the top 10 non-trivial eigenvectors (normalized to unity) {v^ , v[[) 
for each Tf. 

Our aim is to identify an eigenvector for each t/ that can be used to obtain an AIS whose 
corresponding ACSs give a braiding structure similar to the reference case. We denote this 
eigenvector by vl f . By definition, for Tf = 1 the eigenvector that gives the relevant braiding 
structure is the second eigenvector, thus v\ = v\. We calculate the inner product of v\ with 
each of the 10 eigenvectors of the next lowest Tf value, i.e., Tf = 1 — At/, and we compute 
the absolute value of the inner product, /3j = \(vl,Vj Ar/ )|, for j = 2, 11. We select the 
eigenvector Vj Ar/ that gives the highest value of /3j , and refer to this eigenvector as vl Ar/ . 
We denote this value as /3 r/ , i.e., {3 T = \(vl f ,vl f Ar/ )|. Similarly, this procedure is carried 
out for Tf — 1 — 2 At/ by taking the norm of the inner products with vl Ar/ , and so on. 



Using this procedure, we make some observations based on figure [10} The AIS structure 
consisting of 3 ACSs seems to persist until t/ & 0.96. We also observe that as the value of 
Tf decreases further, the structure consisting of 3 ACSs starts breaking up, and at t/ ~ 0.95 
a structure consisting of 16 ACSs can be seen, as shown in figure fl0](c). Further decreasing 



Tf produces another change in the AIS structure, and figure 10(d) shows an AIS with 13 



ACSs. Moving further to the left on the t/ axis, we see AIS structures with 10 and 8 ACSs 



in figures [101(e) and 10(f), respectively. 

This break-up of the AIS structure can be better understood by plotting the f3 T values, 
shown in figure 11, The f3 T value measures the closeness of the eigenvectors along the v* 



branch for two nearby t/ values, and hence captures changes in the eigenvector morphol- 
ogy. We observe that in regions where we can clearly identify the different AIS structures 
consisting of 3, 16, 13, 10 or 8 ACSs, the value of f3 Tf remains almost constant. On the 
other hand, in the transition regions where there is no clearly distinguishable AIS structure 
there is a sharp drop in the f3 T value, signifying a significant change in the structure of 
the eigenvector. Intuitively, this behavior of f3 Tf means that the AIS structure persists for 
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FIG. 11. j3 Tf for a range of Tf values (black curve). The regions where this value stays almost constant are 
the regions of persistence of different braiding structures, separated by transition regions where the value 
drops significantly. The topological entropy of the flow, hf, computed by line stretching, reproduced from 
figure [3j is shown in red. The lower bounds computed by the TNCT for the braids generated by the ACSs 
in that parameter range are shown in blue in the relevant regions. 

a range of Tf and then undergoes a transition (or bifurcation) into another AIS structure. 
Comparing figure 10 with figure 11, we see that the dips in the graph of f3 Tf versus Tf seem 



to accurately capture the boundaries of different AIS structures and transition regions. 

The ACS found using the above method can each be identified with a braid, with the 
number of strands in the braid corresponding to the number of ACS in the domain. For 
parameter ranges over which f3 T remains essentially constant, we find that the number 
of strands in each braid remains constant and the braids are isotopic to each other. We 
also find that dips in the value of f3 Tf correspond to changes in the number of strands in 
the representative braid and changes in the braid topology. We determine the space-time 
braid structure produced by the ACSs in each case by computing the time-shifted transfer 
operators as was done for the braid on 3 strands in the previous subsection. 



Consider the case Tf ~ 0.93, for which the ACSs form a 13-stranded braid. In figure 12 
we show the time-shifted eigenvectors at several different times during the first half-period 
of the flow, and we illustrate the motion of the 13 strands during one full period of the flow 



in figure [L3| The overall structure of the braid consists of four strands each on the left and 
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right of the domain, and five strands in the middle. During the first half-period, strands Bl 
through B4 move together to the right while CI through C4 move to the left, similar to the 
(72 motion of the middle and right strands in the reference case, but strand B5 remains in 
the middle. Similarly, during the second half-period, strand CI remains in the middle, while 
C2 through C4 and B5 move to the left. This motion of the ACS gives rise to the physical 
braid representation shown in figure [l4|(a), which can be written as o^or, where 



Gr = CT_3 G-2 G-3 G-i G-2 G-3 Gq G%Giq GjGqGh Gq Gg G\o G\2 G^G^G\\ (7g G 10 Gq (13a) 

is the braid over the first half period of the flow, and 

GL = GlQ Gn GlO G\2 G\\ G\Q G-4 G _5 G-3 G-Q G-4 G-2 G-7 G _5 G-3 G-i G-Q G-4 G-2 G _5 G-3 G-4 

(13b) 

is the braid over the second half period. All strands return to their original positions after 13 
periods of the flow. This motion is hence topologically distinct from the 3-stranded reference 
braid shown in figures [j](c) and|9](g). The lower bound on the topological entropy of this 
braid is /itn,i3 ~ 0.956, which is lower than /itn,3 for the reference case. 





FIG. 12. AIS identified by second eigenvectors of time-shifted Rt uTf for the first half of the flow period r/, 
where ry w 0.93 (enhanced online). Left to right: Top-row: ti/rf = 0, 0.1, 0.2. Bottom row: ti/rf = 0.3, 0.4, 
0.5. The braiding motion of the 13 ACS is evident, and is illustrated in figure [9] 

Over the range of parameters considered here, 0.84 < tj < 1.0, we find that decreasing Tf 
from unity gives ACSs that generate braids on 3, 16, 13, 10, and 8 strands, successively. The 



physical braid representations are shown in figure [T4j and the corresponding braid words 
are given in the appendix. The basic behavior in all of the bifurcations tends to follow 
the pattern discussed in the 13-strand case — one of the strands from the middle of the 
domain breaks away from the rest of the strands, and hence we get topologically different 
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(a)Initial state 





(b) Movement during pulse 1 
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(c) Movement during pulse 2 



(d) State after 1 period (i.e., 2 pulses) 



FIG. 13. Movement of the 13 ACSs identified from figure 12 over one period of the flow. The ACSs return to 
their initial positions after 13 periods. Note that what was formerly the central set (strand) in the reference 
case now consists of five disjoint sets (strands), shown in red. The strand B5 is 'left behind' in the center 
of the domain during the first half-period and subsequently moves to the left during the second half-period, 
generating a topologically distinct braid from that shown in figure [9] We find this braiding motion of the 
ACS to persists for roughly 0.910 < r f < 0.935. 

braiding in each case. Table [T] lists the topological entropy of these braids, and these values 
are plotted in figure [TTJ We see that the lower bounds provided by the TNCT for these 
braids, /&tn,jvj are indeed sharp lower bounds on the actual topological entropy for the flow, 
hf. Furthermore, /&tn,jv is a good estimate of hf in these cases. This is a key result and 
suggests that ACS-based topological analysis can be used in the framework of the TNCT. 
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Tf range 


AT 




0.960 < T f < 1 


3 


0.962 


0.941 < r f < 0.949 


16 


0.961 


0.910 < T f < 0.935 


13 


0.956 


0.861 < T f < 0.890 


10 


0.936 


0.840 < T f < 0.852 


8 


0.894 



TABLE I. Topological entropy values, /itn,tv, for the braids on N strands generated by the ACSs for the 
various Tf values considered here. 
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(a) (b) 




(c) (d) 

FIG. 14. Physical representations of different braids. The time is increasing from bottom to top. (a) The 
braid on 16 strands, 0.941 < r f < 0.949. (b) The braid on 13 strands, 0.910 < r f < 0.935. (c) The braid on 
10 strands, 0.861 < r f < 0.890. (d) The braid on 8 strands, 0.840 < r f < 0.852. 
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IV. CONCLUSIONS AND FUTURE DIRECTIONS 



We give numerical evidence that the almost-invariant set theory can be applied in con- 
junction with the Thurston-Nielsen classification theorem to find estimates of complexity 
in a wide class of systems. We have used almost-invariant sets as central objects in the 
application of the TNCT by looking at non-trivial braids corresponding to the space-time 
trajectories of these sets. We demonstrate that this procedure can predict an accurate lower 
bound on the topological entropy for a flow even in those cases when the system has been 
perturbed far from one that contains braiding periodic orbits. This work shows that predic- 
tions regarding chaos can be made by considering regions in phase space that move together 
for a finite amount of time. This application broadens the scope of problems where such 
a lower bound can be established and suggests that they should be of use in time varying 
systems that do not have any well-defined periodic points, such as free surface flows, or in 
time-independent systems operating in a parameter regime where there do not exist any 
easily-identifiable fixed points. 

We stress the fact that while the TNCT can be (in theory) applied to any orientation- 
preserving flow (on a two-dimensional domain) that results in a diffeomorphism, the difficulty 
lies in finding the appropriate isotopy class. In the cases where low-order periodic orbits exist 
(and can be found), the isotopy class is defined relative to the periodic orbits. In these cases, 
the stretching produced by the TN representative is a strict lower bound on the stretching 
for any flow in the isotopy class. For the cases in which no obvious, low-order periodic orbits 
exist and the braiding is instead identified via almost-cyclic sets, the predictions of stretching 
are only approximate, since the ACSs do not rigorously qualify as punctures in phase space. 
However, this ACS approach identifies a small set of 'approximate punctures' that give an 
accurate representation of the topological entropy for the cases we have considered. 

The class of set-oriented transfer operator methods discussed in this paper have been 
applied to systems accessible solely by time series dataP^, systems with stochasticity 62 , and 
non-autonomous systems defined by data over finite time, where finite-time coherent sets be- 
come the objects of interest 63 . Hence, the connection between set-oriented statistical methods 
and topological methods in dynamical systems made in this paper provides an additional 
tool for analyzing complex systems, including those defined by data. Ideally, one would want 
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to apply the notion of braiding coherent sets to an experimental setting, where one could 
ask questions regarding the role of braiding coherent sets in complex physical systems. For 
example, what implications does large-scale braiding have for stirring in geophysical flows? 

There is work on spectral analysis of the "mixing matrix"^ that is closely related to the 
identification of almost invariant sets and ACSs, and connections have been made between 
this mixing matrix analysis and the "strange eigenmode" that arises from spectral analysis 
of the continuous advection-diffusion operato]P^. The topological information available 
from analyzing trajectories of ACS suggests that a similar approach can be applied to the 
braiding of eigenvectors in these related methods. 

We note that for the lid-driven cavity system, some dynamical 'memory' of the saddle 
point stable and unstable manifolds remains even where there are no more saddle points] 
compare figure [6| where such objects exist, with the wisp-like features of figures [7] and [8j 
where they do not technically exist. Transfer operator methods seem to reveal structures 
resembling stable and unstable manifolds of invariant sets, but these are now associated with 
almost-invariant sets. Given the importance of stable and unstable manifolds in applications, 
future work could extend this concept of stable and unstable 'almost-invariant manifolds' 
to systems that are not exactly periodic, and possibly far from periodic. For example, in 
applications to stability and controP^23, 'almost-invariant manifolds' may prove useful for 
geometric interpretations of stable and unstable directions in the phase space. They may also 
help in making connections with invariant manifold-like objects detected by other methods; 



cf. figure 15 



Another intriguing possibility is opened up by considering the spectral dependence of 
transfer operators. Even in relatively simple systems such as shown in figure [lOj one sees 
interesting changes in the eigenvectors of the discretized Perron-Frobenius operator as a 
system parameter is varied. Some modes increase in importance while other modes decrease. 
Different modes can correspond to dramatically different behavior. This observation may 
have an interesting application to systems defined from data: prediction of dramatic changes 
in system behavior based on mode variations. That is, a system may contain hints of 
a bifurcation before a bifurcation occurs, and this premonition may be teased out via a 



transfer operator approach, similar to 59 . To take a vivid example, consider figure 15, where 



we show the splitting of the ozone hole in 2002 71 , which we might classify as a 'Duffing 
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bifurcation', where the bifurcation parameter is some physical quantity changing slowly 
with time. Given only the weather observations up to some time, £ now , before the split, e.g., 
September 20, 2002 (figure [l5](a,b)), might the spectrum of the discretized Perron- Frobenius 
operator reveal an emerging, and topologically distinct, eigenvector which is increasing in 
importance, soon to become the dominant mode (figure [l5|(e))? To determine trends in 
the change of spectrum of a transfer operator in real-time, one may need to consider the 
possibility of using the infinitesimal generator of the discretized Perron-Frobenius operator, 
which could give several orders of magnitude increase in speed for determining the spectral 
characteristic^. 
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FIG. 15. (a,c) Atmospheric Lagrangian coherent set boundaries (LCSBs) for days surrounding an Antarc- 
tic polar vortex splitting event in September 2002 (based on NCEP/NCAR reanalysis data); attracting 
(repelling) curves are shown in blue (red). Before and after the splitting event in late September, we see 
an isolated blob of air, bounded by LCSB curves, slowly rotating over Antarctica. The vortex pinches off, 
sending the northwestern part of the ozone hole off into the midlatitudes while the southwestern portion 
goes back to its regular position over Antarctica. We note that the time-varying structures that bound the 
ozone hole fragments provide a framework for understanding the geometry of atmospheric transport; they 
act as analogs of stable and unstable manifolds of saddle-like almost cyclic sets. (b,d) The corresponding 
daily ozone concentration (based on NASA TOMS satellite data). (Adapted frorrl^) (e) Shown in schematic 
is a possible change in ordering of the leading eigenvectors for the polar vortex, suggesting a new means 
to determine global bifurcations of a system. Even when the 'single-blob' eigenvector is dominant, there 
may be 'double-blob' eigenvector whose eigenvalue is increasing. Somewhere near the crossover point, the 

'double-blob' eigenvector will be dominant and determine the large-scale characteristics of the flow. 
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Appendix A: Braid words for different braids 

We give braid words for two pulses (i.e., one complete flow period) in each case. Here, a R 
refers to first half of the time period and a L refers to second half. The complete braid word 
for each braid is formed by composing the two words, i.e., o^or. 

• Braid on 3 strands. 
ctl = cr-i 

• Braid on 16 strands. 

= ^12^13^12^14^13^12^15^14^13^12^-5^-6^-4^-7^-5^-3^-8^-6^-4^-2 

a_9a_7a_5a_3a_ia_ 8 a_ 6 a_4cr_2cr_7a_5a_3a_ 6 a_4a_5 

• Braid on 13 strands. 

• Braid on 10 strands. 

a L = a_3Cr_4Cr_2Cr_5Cr_3Cr_ia_4a_2Cr_3Cr_iCr_2Cr_i 
a R = (Jg(J%(Jg(J 7 (J%(jQ(Jg(J 7 (J5(J%(jQ(J 7 

• Braid on 8 strands. 

(j L = a_ia_3a_4a_2cr_ 3 cr_ia_2 
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